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We study mean field dynamo action in renovating flows with finite and non zero correlation time 
(r) in the presence of shear. Previous results obtained when shear was absent are generalized to 
' the case with shear. The question of whether the mean magnetic field can grow in the presence 

of shear and non helical turbulence, as seen in numerical simulations, is examined. We show in a 
, general manner that, if the motions are strictly non helical, then such mean field dynamo action is 

^Nj ' not possible. This result is not limited to low (fluid or magnetic) Reynolds numbers nor does it use 

> I ' any closure approximation; it only assumes that the flow renovates itself after each time interval r. 

^ |. Specifying to a particular form of the renovating flow with hehcity, we recover the standard dispersion 

■ relation of the c?Q. dynamo, in the smaU r or large wavelength hmit. Thus mean fields grow even 

' in the presence of rapidly growing fluctuations, surprisingly, in a manner predicted by the standard 

, quasilinear closure, even though such a closure is not strictly justifled. Our work also suggests the 

possibility of obtaining mean fleld dynamo growth in the presence of helicity fluctuations, although 
having a coherent helicity will be more efficient. 



> 



PACS numbers: 47.27. W-,47.65.Md,52.30.Cv,95.30.Qd 
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' I. INTRODUCTION 



Understanding the origin of coherent large scale magnetic fields, observed in astrophysical systems from stars to 
galaxies, is of fundamental importance in astrophysics. The standard paradigm invokes the amplification of small seed 
magnetic fields due to dynamo action, typically involving large scale shear flows combined with helical turbulence 
^ ^ [IIj HI- Recent numerical simulations have also raised the possibility that coherent fields can arise in non helical 
] turbulent flows in the presence of shear 0] , although the mechanism of how this happens is unclear 0, [1] . 

The dynamics of large-scale magnetic fields is generally described by using the equations of mean field electrody- 
namics. Here one defines the mean magnetic field B and mean velocity field U by suitable averaging over the small 
. scales corresponding to the turbulent fluctuations. The mean field evolution is governed by the averaged version of 
the induction equation, which results in an extra term, the turbulent electromotive force £ = {u x b), where u and 
b are the fluctuating velocity and magnetic fields. Expressing £ in terms of the mean fields themselves is a closure 
problem, even for prescribed velocity fields. If the fluctuations can be assumed to be small, one can employ the 
quasi-linear approximation or what is traditionally known as the first order smoothing approximation (FOSA) to 
5^ ■ express the turbulent electromotive force £ in terms of a term proportional to B (the a~effect) and one proportional 
to the mean current density (turbulent diffusion). The a-effect, which depends on the helicity of turbulence is crucial 
to amplify mean fields. 

However, turbulent motions above some modest magnetic Reynolds number lead to a fiuctuation dynamo and a 
rapid growth of magnetic noise P, 0i ■ The growth rate of the fluctuation dynamo is typically larger than the growth 
rates associated with the mean-field dynamo. In the presence of this rapidly growing magnetic noise, the validity of 
the quasi-linear approximation or FOSA becomes suspect. The question then arises whether one can indeed make 
sense of mean field concepts like the alpha effect. 

In this context, considering exactly solvable fiow models becomes very useful. For example, if one assumes a random 
flow whose correlation time is exactly zero, dynamo action, on both small and large scales can be studied in great 
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analytical detail [110; 13; however, such a flow is unphysical. Renovating flows discussed by several authors (and 
references therein) provide, on the other hand, models involving flows with a finite non zero correlation time but which 
are still analytically tractable. In renovating flows time is split into successive intervals of length r and the stochastic 
component of the velocity in the different intervals are assumed to be statistically independent realizations of an 
underlying probability distribution (PDF). As the flow loses memory between different time intervals, the evolution of 
the moments of the magnetic held over any one time interval can be calculated by averaging over the underlying PDF. 
Considering random helical renovating flows, Gilbert and Bayly (GB) [l^] showed that the magnetic fleld becomes 
increasingly intermittent with time. Nevertheless, the mean magnetic fleld can still grow with a growth rate which 
approaches that of a standard mean-fleld o? dynamo in the limit of small r [1, [l^l . These works thus provide explicit 
demonstration that the growth of magnetic noise need not destroy the growth of the mean fleld even in the case of 
flows (which have this periodic loss of memory) with flnite non zero correlation times. 

In the present work we generalize some of these results to renovating flows incorporating also a large scale shear. 
Our primary motivation is to examine if the introduction of shear can lead to a mean fleld growth even if the stochastic 
velocity fleld is non helical. It turns out that, once properly formulated in terms of renova ting shearing waves, the 
details of our calculation share some crucial features with those of Gilbert and Bayly (GB) [101 ■ However, GB have 
not given these steps; so to aid the presentation of our results, we begin in § 2 with a presentation of the main results 
of [lO| for the mean fleld evolution in renovating helical flows. In § 3 wc formulate the problem of renovating flows 
with a background linear shear, and prove a general result that there is no dynamo action when the flow is strictly 
non helical. We also consider a particular example of renovating shearing waves with helicity, caused by overdamped 
external forcing and recover the dispersion relation for the a^f2 dynamo. The flnal section summarizes and presents 
a discussion of our results. 



II. RENOVATING HELICAL FLOWS WITHOUT SHEAR 

We re~derive here the results of Gilbert and Bayly on mean fleld evolution in a model helical renovating flow, 
in the absence of background shear. Consider the induction equation for the evolution of the magnetic fleld, 

— = V X [u X B - r;V X B] . (1) 

GB assumed the velocity fleld u to have zero mean with only a turbulent component. In each renovating time interval 
T, they took 

u(x) = a sin(q • x + ?/') + cos(q • x + -0) . (2) 

with the conditions 

a • q = 0, b = qxa/g (3) 

which implies incompressibility, that is V ■ u = 0. The parameter h satisflcs — 1 < /i < 1 and determines the helicity 
of the flow. This helical flow is made random by choosing the parameters of the flow randomly and independently 
from an underlying PDF, for every time interval r. The ensemble considered is the following: In each time interval 
\(n — 1)t, nr], (i) -0 is chosen uniformly random between to 27r ; (ii) the propagation vector q is uniformly distributed 
on a sphere of radius q ; (iii) for every flxed q, the direction of a is uniformly distributed in a circle of radius a in 
the plane perpendicular to q. The parameters (g, a, /i, r) are non-random and completely describe the renovating 
flow. The randomness of -0 in condition (i) ensures statistical homogeneity, whereas conditions (ii) and (iii) ensure 
statistical isotropy of the flow. 

The evolution of the magnetic fleld from time (n — l)r to ut is given by 

B,(x, nr) = y g„ (x, y)B, (y, (n - 1)t) d^y (4) 

where Qij (x, y) is the Green's function of the induction equation Eq. ([IJ . Qij is random due to the randomness of the 
turbulent velocity fleld. We take the ensemble average of this equation over the ensemble described above and note 
that the velocity in any time interval \(n — 1)t, nr] is uncorrclatcd with the initial magnetic fleld at time (n — l)r. Thus 
the average of the product of Qij and Bj can be written as the product of the averages, an important simpliflcation 
arising from the loss of memory of renovating flows. The mean fleld B then evolves as, 

B,(x,nr) = [g~i^,y)B,iy,in-l)T)d^y. (5) 



3 



Further, from the statistical homogeneity of the renovating flow, one has Qij{x,y) ~ Gij{x — y); then Eq. ((5)), which 
is a convolution in physical space, becomes a product in Fourier space. Defining the spatial Fourier transform. 



B,{k,t) ^ J d^x B,(yi,t)e~'^'' (6) 

we have in Fourier space, 

t,(k, nr) = Gy(k) X tj{k, {n - 1)t) (7) 
where the response tensor, G'ij(k), is defined by, 

G,,(k) = Jg~{x-y)e-'''-^^'yWx. (8) 

The mean magnetic field will grow exponentially if its Fourier component is a eigenvector of the matrix Gij (k) with 
eigenvalue cr, whose magnitude is greater than one. For such an eigenvector we have 

ii (k, nr) = a" B, (k, 0) or Bi (k, t) = cr*/"" B, (k, 0) (9) 
and so the growth rate of this eigenmode is given by 

A=i^. (10) 

T 

Hence, the response tensor G'ij(k) contains all the information about the growth or decay of the mean magnetic field. 
We now proceed to calculate it for the renovating velocity field of Eq. 



A. The response tensor 



In order to explicitly calculate the response tensor for the evolution of the mean field in the renovating helical flow, 
an important simplification was introduced by GB. The renovation time r was split into two equal sub-intervals. In 
the first sub- interval (step 1) resistivity was neglected and the field was just frozen-in and adverted with the fiuid, 
with twice the original velocity. In the second sub-interval (step 2), induction by the velocity was neglected and the 
field was assumed to diffuse with twice the resistivity. Although such an assumption seems plausible in the limit of a 
short renovation time — as this would be one way of numerically integrating the induction equation — GB did not give 
any rigorous justification. For the present purpose, we adopt the same simplification as GB. Thus we consider the 
evolution of the magnetic field in these two steps and then Fourier transform the resulting averaged Green function. 
Step 1: During the time interval to r/2 we assume rj = 0, and double the value of the velocity field. Then Eq.([T]) 
becomes just the ideal induction equation, 

with the standard Cauchy solution given by 

B^r.t) = -^B,(y,to) = Jy(r-)B,(y,io). (12) 

Here B(y,io) is the initial magnetic field; r{t) is the position of a fluid element at time t, which was originally at a 
'Lagrangian' position y at time <o- Note that the fluid elements follow the integral curves of the velocity field, with 

dx/di = 2u = 2asin($) + 2b/icos($) (13) 

where we have substituted from Eq. ([2]), assumed twice the velocity for step 1 and defined the phase 4> = q • x + -0. 
From the incompressibility condition, we have (i$ / dt ~ 2(\ - vl — Q and thus Eq. (fT3)) can be integrated to give at time 

t = r/2, 



r 



= y + ru 



— y + ar sin (q • y + V') + brh cos (q • y + ^) ■ 



(14) 



4 



Here we have used the constancy of the phase 4> and set it equal to its initial value $ = q • y + -0. Thus the Jacobian 
is 

Jij{r) = TT^ = fli^jT cos (q • y + V>) - &j(7j7ir sin (q ■ y + t/") (15) 

Step 2: During the time interval t/2 to t the turbulent velocity field is zero and there is only diffusion present, with 
a resistivity 2?7. The induction equation then reduces to a diffusion equation for the magnetic field as 

^ = 2v^'B (16) 

The solution of this equation is given in terms of the resistive Green's function 

(x-r) 



4rjT 



(17) 



The total Green's function defined in Eq. Q is simply the product of the two Green's functions in the above two 
steps 

a.,(x,y) = g"(x-r(y))J„-(r(y)), (18) 
where we have written explicitly r as a function of y. The response tensor defined in Eq. ([8]) then becomes 



where in the second step we have done the integral over x. The overhead bars in Eq. (|19p denote ensemble averages, 
and we will see below that due to the statistical homogeneity of the renovating flow, this averaged quantity does not 
depend on y, but only on k. Also note that in typical astrophysical systems, the resistive timescale will be much 
larger than the renovation time and the value of ?yfc^r is typically much smaller than unity, and so can safely be set 
to zero. A non zero but small ry will decrease the growth rate by a negligible amount. 

We now evaluate the ensemble average in Eq. (fT9|) . Note that GB state the final result, omitting all intermediate 
steps. We give the detailed steps in Appendix |X] since they are of use in the case when shear is present. Here we list 
the important intermediate steps and the final expression. Let the angle between k and q be 9; we will treat this as 
a colatitude and denote the azimuthal angle of q by (j). Let the component of k perpendicular to q make an angle (jj 
with a (see Fig. [!}. On averaging over phase we get 



Gijik) = Sij Joirakx sin 9) —— [e„nnhnqnqj] Jiirakx sin 9) (20) 

Xksmu 

where x = (cos^ 4>+h'^ siv? 4>)^l'^ and the overhead bars denote ensemble averages over the remaining random variables 
0, (/) and 0. Averaging over direction of q (i.e. averaging over the angles 9 and (/), we have 

Gij(k) = Sij go{Tak,h) + ^^'^"''^^'^"'^'^ gi{Tak, h) , (21) 

where 



sm(sx , ,^ 1 ( sm SX cos(sx)\ 
go{s,h) = ; gi[s,h) = - - — (22) 

SX X V (sx) SX / 

and now the overhead bars denote ensemble averages over the random variable (f> (for maximally helical flow with 
/i = ±l,sox=l, and the response tensor becomes independent of the random variable (f>. In the rest of this paper 
h can take any value between —1 and +1). We recall that in step 1, we evolved the field for a time interval t/2, but 
with twice the velocity (i.e. 2a instead of a); nevertheless the combination ar which appears in the above equation 
remains unaffected. 

One of the eigenvectors of Gij is k with eigenvalue unity. But this eigenvector can be ignored since the magnetic 
field mode must be orthogonal to k (i.e. we must have V • B = 0). The relevant eigenvectors of Gij are found to be 
(— z, 1,0)"^ and (z, 1,0)'^, with the corresponding eigenvalues, (j+ and (t_ , given by 

raqh 

<7± ^ go[Tak,h) =F -^^gi(jak,h) (23) 



FIG. 1: The difTerent angles defined in tiie text. Tiie angle between k and q is 9. The vectors a, b and q are mutually 
perpendicular to each other. The component of k perpendicular to q makes an angle (j) with a. 



The growth rates A±, are given by 



Ah 



ln(a±) 



In 



goiTak,h) =F ^^^^gi(Tafc, /i) 



(24) 



Here we have divided by the full time interval r to get the growth rate. Following GB it is readily verified that 
both ± mean field modes can grow for sufficiently large renovation times. GB also demonstrated that the magnetic 
field becomes increasingly intermittent, in the sense that higher order single-point moments of the field grow faster. 
Therefore, as advertised, the mean field dynamo operates efficiently in this case even in the presence of strongly 
growing magnetic noise. It is of interest to look at the growth rate in the limit of small renovation times such that 
akr <C 1 (this limit can also be looked upon as a long wavelength, or small k, limit at a fixed r), when 



Tk 



ha?Tq 



(l + fe^)a^ 
12 



±ak — rjtk 



(25) 



where a = — (l/3)u- (V x u)(t/2) and 7^4 = (l/3)u • u(t/2) are the turbulent transport coefficients.^ Thus the growth 
rate for the case of small akr is identical to the growth rate of the standard mean-field dynamo usually obtained 
using the quasilinear approximation or FOSA. We now turn to consider the influence of shear. 



III. RENOVATING SHEAR FLOWS 



We now investigate the evolution of the mean magnetic field in scenarios when there is a mean shear flow over 
and above the background turbulence. Shear flows and turbulence are ubiquitous in astrophysical systems. Recent 



In these expressions for a and r]t, the factor t/2 appears instead of r, because the transport coefficients are really time integrals of 
correlation functions; see Eqs. (6.19) and (6.20) in [y. For renovating ffows in which the velocity field u acts over the full interval r 
(i.e. if we were dealing with Eq. [inductioneqn] , without the two-step prescription of GB), this implies averaging t over the time interval 
(0,t), which is equal to t/2 . However, we get the same result even for the two-step prescription of GB, because two effects contribute 
in precisely opposite ways: when the velocity field is doubled in value, the transport coefficients quadruple because they are quadratic 
in the velocities; however, the doubled velocity field is ON for only the first half of any time interval r. Hence, averaging "i" over the 
time interval (0, t) now implies integrating t over the time interval (0, r/2) and then dividing it by the full interval t, which gives r/8. 
Thus we obtain 4 X r/8 = r/2 . 
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work suggests that the presence of shear may open new pathways to the operation of large-scale dynamos [3|-|5| . For 
SHTiplicity we consider the background mean velocity to be a linear shear flow. Let (ei, 62, 63) be the unit vectors of 
a Cartesian coordinate system in the lab frame, x = {xi,X2,X3) the position vector. Without loss of generality, we 
choose this mean velocity to be in the 62 direction and varying linearly with xi. Thus the velocity field is given by 

U{x,t) = Ush + Uturb = Sxie2 + Uturb (26) 

where S is the constant rate of shear coefficient. The turbulent velocity field Uturb is composed of renovating shearing 
waves with quite general amplitudes at this point. In particular we take 

uturb(x, t) = A(i, q) sin [Q{t) ■ X + ^A] + C(t, q) cos [Q{t) ■ x + V] , (27) 

where the wavevector is a shearing wavevector of the form Q = ((5i,Q2,<93) = (91 — Sq2{t — ti), 92, 93), and q ~ 
(gi,g2,(?3) its initial value at the beginning of each renovation period, i.e at i = t^. Note that q is chosen randomly 
from a specified PDF (sec below) for each renovating period. We will see that such a form of Uturb naturally arises 
when we consider Fourier modes of the velocity which satisfy the momentum equation in a background linear shear 
flow. We will also later adopt explicit forms of A(t, q) and C(t,q); but several of the conclusions that we arrive at 
are quite general insensitive to the functional forms of A and C. We also assume the turbulence to be incompressible 
with V • Uturb = 0, which implies 

Q(t)-A(t) ^ 0; Q(i)-C(t) (28) 

Thus the amplitudes have to shear in an opposite sense to the wavevector so as to maintain incompressibility. The 
shearing wavevector can be written in a compact form as Qj ~ Qi7iji^(t — ti)), where jij(t) is the shearing matrix 
defined by, 

7,,(t) =. % + 5,26jiSt. (29) 
The helicity of the turbulent velocity field Uturb (x, t) is 

i^ = Uturb-(VxUturb) = C-(Qx A), (30) 

and this vanishes unless both C or A are nonzero. 

As in the previous section we consider the turbulent flow to be a pulsed renovating flow. The turbulent velocity 
field is assumed to be ON for a time interval t/2, with twice its amplitude and with diffusion absent. For the next t/2 
interval, the turbulent velocity field is OFF and only the diffusion is present with resistivity 27]. On an average, the 
turbulent velocity field is then correlated only for a time interval r. The mean shear fiow on the other hand is always 
present for the full time interval r. The turbulent flow is randomized by considering an ensemble similar to that 
assumed for the renovating flow without shear. In each time interval [{n — 1)t, nr] (i) is chosen uniformly random 
between to 27r (ii) the propagation vector q is uniformly distributed on a sphere of radius q. The randomness of 
A is decided by the explicit form of A itself, which we fix later, when solving for the explicit form of the response 
tensor. For the general analysis we will not require it. The response tensor will be essential here too to determine the 
growth or decay of the mean magnetic field modes. 



A. The response tensor 



We compute the response tensor for the evolution of the mean magnetic field again in two steps. 
Step 1: During the time interval t = to t = t/2, rj = Q and Eq. ([T]) reduces again to the ideal induction equation, 
whose solution is as before the Cauchy solution 

B,{f,t) = ^Bj{y,Q) = J,,{f)B,{y,0) (31) 

Here f gives the location of the fiuid element at time t, which was at time t = at the location y. These positions 
are now on the integral curve of the sheared and turbulent velocity field, and so the trajectory now obeys 

dx/dt = Ush + 2uturb = Sxie2 + 2A(t,q)sin$ + 2C(t, q)cos$. (32) 

We have substituted here from Eq. (|26p . assumed twice the turbulent velocity for step 1 and defined the phase 
$ = [Q{t) • x + -0]. Note that we have not doubled the shear velocity, as we keep the shear flow throughout the full 
period (0,t). From the incompressibility condition, we have Q • Uturb = 0. Therefore 



d^/dt Q ■ X + Q(i) ■ X = -S'g2a;i + Sq2Xi + 2Q • Uturb = 0. 



(33) 
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The constancy of $ can be used to express it in terms of the initial position of the fluid element y, and the initial 
wavevector Q(t = 0) = q, that is we can write $ = q • y + -0. Then Eq. ([5^ can be integrated to give 



Tj = aj(t, (7fc)sin(q • y + V) + Cj(i,(7fe)cos (q • y + + ?/j 
where Tj is a sheared position vector that will be of use later, and the coefficients aj and Cj are defined by, 



(34) 



hj{t,q,) = ljp{~t) / 2Ap{t',q,)dt' + SSj2 



t pt' 



2Ai{t",qi)dt"dt' 



Jo 
t rt' 



Cj{t,qi) = 7jp(-t) / 2Cp{t',q,)dt' + SSj2 f f 2Ci{t" ,q,)dt"dt' 
Jo JO Jo 



Therefore the Jacobian matrix is, 

dfi 



J^j{r{_y),t) = — = "fip(t)[Spj + d.p{t)qj cos {q ■ y + ip) - Cpit)qj sin {q ■ y + tp)] . 



(35) 



(36) 



Wc need this Jacobian to be evaluated at the time t = t/2. 

Step 2: During t/2 to t the turbulent velocity field is zero and there is diffusion present along with shear. The 
induction equation then reduces to the following form 



d_ 

dt 

The sheared Green's function for this equation is Q 

7., (t/2) 



u,h • V B = (B • V)u,h + 2r]V^B 



(det|D| 



a/2 



exp 



4:T]T 



(37) 



(38) 



where Xi = 7y(— r/2)a;j — f ;, and Dij is a symmetric matrix whose inverse is given by 

1 -St/4: ■ 



-St/4: 1 + (5t)2/12 
1 



(39) 



Note that in the above computation wc have taken the initial time to be i = and the final time as t = r; but 
the same steps and calculations go through for any renovating time interval [{n — 1)t, nr], with the initial time being 
t ~ (n — 1)t, the final time being t = tit, and the initial wavevector q = Q((n — I)''')- Hence during any such time 
interval r, the magnetic field at time t = ut is related to the magnetic field at time i = (n — 1)t by 



S.(x,nr) = j gl(^^,fj-^ j,^. (^(y),0 (y, (n - l)r) 



(40) 



Wc would like to calculate the response tensor starting from the above evolution equation. For this we first define the 
Fourier transform of Bi(x, t) by expressing it in terms of the shearing waves, 



S,(k,t)- j S,(x,i)e-'^(*) ''d3 



(41) 



where we have defined the sheared wavevector Kj{t) = kijij{—{t — ti)) = (fci — Sk2{t — ti), k2, k^) and k = K(ii) is 
the initial wavevector at time t^, which for each step we take to be the time (n — 1)t. We will see that the evolution 
of the mean magnetic field is especially transparent when the field is expanded in terms of shearing waves in Fourier 
space. 
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We take the Fourier transform of Eq. (|40)) . and change the integration variable from f to y. The Jacobian for this 
transformation is unity as the flow velocity which maps y to r is divergence free i.e. (Pf = (Py. We get 

B,(k,nT) = y"B,(x,nT)e-''^(""'" d^x 



(42) 



where we have also expressed the initial field at time {n — 1)t in terms of its Fomier transform. To do the integral 
over X, we use the identity ^jm {t)^mp{-t) = 6jp to write 

exp (-iK • x) = CXp {-iKj^jr,iiT/2)[jrnpi-T/2)Xp - rrn] ^ iK j^.jyn{T /2)f„i]) 

= exp [-i/Cj7_,-™(r/2) (x™ + r™)] (43) 

and change the integration variable from Xm to x„i ~ jmp{~T/2)xp — r^. Since detj^j ~ 1, we can write dPx — dPx. 
The integration over x now becomes a Fourier transform of the sheared resistive Greens function Q and we get 

/ Ql (x, r, e'^^-d^x = J 0^ (detlDI) exp (^^^^^) e---/^-e-- K/^^]d3f 

= 7^fc {r/2) exp {-r?T/C„7„„(T/2)7^-i7,p(T/2)A-,} e^^^'t^-z^'^] 

= GMe-'^'' (44) 
where = Kj^jm{T /2) = fcj7j>„(— t/2). Then we have 

i3,(k, nr) = 7,fc (r/2) exp (-i^rlpD^yK) ( ( J,., (f (y), ^) ^ (1, (n - l)r)e'' y e-^^'^ . (45) 



(27r) 

Taking the ensemble average of Eq. (|45p , we get for the mean field evolution, 



i?,(k, nr) = 7,fc (r/2) exp (-Tyrfepi^-ifc,) ^ ^ e-'i'-Jfe, (f(y), ^)e*'yd3y 



S,(l,(n-l)r)^ (46) 



Here we have assumed as before that the velocity field during the time interval (nr, (n — l)r) and the initial magnetic 
field at time (n — l)r are statistically independent. Note the dependence on the stochastic parameters comes only 
through f{y). As fcm = fcj7j>„(— r/2) and r,n = ImjiT /2)rj, we have k • f = k • r and the quantity in the square 
bracket can be written as 

Je-^^ -Jkj (f{y)J-^e'' ySy = j e-^^ i-v) Jk, (r(y), ^) e-^^ ye^' ^d^y = e-'k.(r-„) j^^. (r(y), (27r)3<5(l - k). 

(47) 

In the last step we have used the fact that the averaged quantity is independent of yi , as can be easily seen by doing 
the averaging first over the random phase of the turbulent field; this also follows from the statistical homogeneity 
of the turbulence. The yi independence is explicitly shown by calculation below. We then have for the evolution of 
the mean field, 

nr) = Gy (k);^j (k, (n - l)r) (48) 

where the response tensor dj (k) is now 



G'ij(k) = 7,;; (r) + a;5j cos(q • y + i/O - c/g^ sin(q • y + -(/')] e ''^'^'^ y) exp (-jyrfcpAfp^fcr) (49) 

where Mij = 7jm(— r/2)D~p7jp(— r/2) and the coefficients a; and c; are shorthand for aj(r/2, q) and c;(r/2, q), given 
in Eg. pSI) . One can see that form of the response tensor in Eq. (j49p with shear is similar to the form of response tensor 
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in Ea. (|T9|) without shear and reduces to the latter when S = 0. Similar to the case without shear, we see that the 
effects of resistive dissipation appear only as a separate exponential term. Since it is small in astrophysical systems 
of interest, henceforth we set 77 = 0. 

We now average over the phase ip of the turbulent velocity field, which leads us to the following expression 



cri] (k) 



5, Jo (V(k.a)^ + (k.c)^) - .^^#=4^fi^ Jr ( V(k • a)^ + (k • c)^ 
^ ^ ■^/(k • a)^ + (k • V 



(50) 



where the overhead bar now refers to the averaging over the directions of q and over the randomness of A. One can 
reach some conclusions about the decay of the mean magnetic field at this stage of the averaging itself. The mean 
magnetic field evolves as 



B,{k,nT) = 7,i(T)t7,,(k)Bj(k,0) (51) 

The growth or decay of the mean field mode is governed by the product of the two matrices, the shearing matrix 7(t) 
and the shear-turbulence tensor <T(k). It is known that the field grows linearly due to the continuous shearing of the 
background fluid which causes the i3i(k) component of the field to be continuously advected along the 62 direction. 
This is reflected by the presence of the shearing matrix 7(r) in the expression of the response tensor G(k) and is 
completely natural as well as expected. The turbulent stretching of the field lines due to the transfer of energy from 
the turbulent pulses of the fiuid can too lead to the growth of the mean field and the shear-turbulence tensor cr(k) 
contains precisely this information through its dependence on the random parameters. Hence, it is of interest to look 
at the structure of the cr(k) tensor and calculate its eigenvalues depending on which the mean field grows or decays 
exponentially. 

When helicity is zero, then either a or c is zero from Eq. ((30|) and Eg. ((35)) . Then the second term in cr/j(k) vanishes 
and we get 



(Ty(k) = %Jo(k-a) (52) 



where we have chosen c = without any loss of generality. In this case, the eigenvalue of <T(k) is just a = Jo {k ■ a) . 
Note that the maximum value of the Bcssel function Jq (x) is unity. Hence after averaging over all the possible values 
of k • d as per the ensemble chosen, we must necessarily obtain Jq (k ■ a) < 1. This shows that in the absence of 
helicity the mean field modes eventually decay with a decay rate A = (T)~"'^logcr (see Eg. pO]) '). Therefore, quite 
generally, there is no mean field dynamo if the turbulent velocity is strictly non helical, even in the presence of shear. 



B. Forced overdamped shearing wave 

We now solve for the form of the <Jij(k) by taking a particular form of the Uturb in Ea. (|27)) obeying the following 
forced, damped Euler equation: 



d_ 

di 



Sxi 



_d_ 

dx2 



Uturb + S'ltJui-bez + (uturb ' V)utuib = -Vp - 



Uturb 



(53) 



where is a given damping time and f(x, t) is the external forcing which is assumed to satisfy V • f = 0. In the 
approximation dut^Th/dt ^ Uturb/'''d, the wave is assumed to be overdamped, saturating quickly in time Td to its 
terminal velocity. In Ea. (|B13[ ) of Appendix [B] the following solution is derived: 

Uturb(x, t) = A(t, q) sin [Q(t) ■ x + tA] + C{t, q) cos [Q(t) • x + ■i/'] , where 



Ci,: 



ai,3 + STdai 



- STdQiQ2 



a2 + Stai — STdai 



Ql + Ql 



h { C1.3 + STdCi 



Qi.iQi 



- STdQlQ2 



C2 = h{c2 + StCi - STdCi 



Q^-STdQiQ2_ 

Ql + Ql 



- STdQiQ2 



(54) 



where Qj — qi'-fij{—t), q-a = 0, c = qxa and h such that — 1 < ^ < 1 determines the helicity of the flow. Here 
the forcing function is related to the constants {01,02,03} and {01,02,03}. In the limit of very strong damping, i.e. 
— ^ 0, the above expressions simplify and can be written compactly as 



A,. 



(55) 
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The turbulent velocity field now becomes 

[uturb(x, t)]^ = j,j {t)aj sin [Q(t) • x + -0] + /i7„- {t)cj cos [Q(t) • x + (56) 
Substituting Eg. ([55)1 in Eg . ([55)) . we find that a = ar and c = her ; hence the response tensor in Eg . ([50)) becomes 



Gy(k) = 7,Kr) 



6, Jo (V(k.ar)^ + (k.cM^) - ^ ^^^^^§i=|=#^ ( V(k • ar)^ + (k • c^^ 
^ ^ ■^/(k • a)^ + (k • c/i)^ ^ 



(57) 

As before, let the angle between k and q be 9\ we will treat this as a colatitudc and denote the azimuthal angle of q 
by (f). Let the component of k perpendicular to q make an angle <j) with a. Then, on averaging over the phase ip we 
can write 



G.y (k) = -fii{T)\6ij Ja{Takxsu\e) - [timnkrnqnqj] Ji (rofcx sin g) I (58) 

where x = (cos^ cft+h'^ sin^ 'PY^'^ ^md the overhead bar denotes ensemble averages over the remaining random variables 
6, (p and (/). Comparing with Eg. (1^01) . we see that the form of the response tensor is identical to its form when shear 
is absent with the only difference being an overall ^uir) factor. Since the ensemble average is done at one arbitrary 
instant, the final form of Gij(k) is identical to Eg. ipTj) with the extra factor. Hence, we have 

G„(k) = laiT)crijik) 

aij{k) = Sij goiTak,h) + ^^'^"■^^rmkm ^^^^^ 



sin(5x) , , s 1 / sin(sx) cos(sx)\ .rn^ 

9,{s,h) = ——■ g,(.,/,) = (59) 

and now the overhead bars denote ensembles average over the random variable </> (for maximally helical flow with 
/i = ±l, sox = l, and the response tensor becomes independent of the random variable (f). The shear turbulence 
tensor a has the two non-trivial eigenvectors (—1,1,0)"'" and (z,l,0)'^, with the corresponding eigenvalues, a+ and 
(7- given by 

(J±^ go{Tak,h) =F ^^^y^5i(raA;, /i) (60) 

For zero helicity, the second term vanishes and dynamo action is absent, as was shown in the general case in the 
previous section. Moreover, even if there were mirror-symmetric fluctuations in /i, this would not lead to a dynamo. 
This is because go{s, h) and 51 (s, h) are even in h, while the co-efficient of the second term of the response tensor in 
Eq. (|59p is linear (and hence odd) in h. Thus on averaging the response tensor over any symmetric PDF of h with 
zero mean only the first term of Gij survives and there is no dynamo. This conclusion is similar to that obtained by 
GB for the case without shear. 



IV. a^Q DYNAMO 



Let us look at the expression in Eg. ([5^ of the response tensor in the case of small correlation times, when akr <^ 1. 
Then to quadratic order in r, we get 



G,,(k) = 5, 



1 - ^^^(rafc)2 



+ Si2SjiST + i ti-jmkm 



6ij{l-'i]tTk ) + 6i2SjiSt - iTae^jmkm 



(61) 



where a = —(1/3) Uturb • (V x Uturb) (^/S) and rjt = (1/3) [u • u] (t/2) are the turbulent transport coefficients. The 
mean field evolution equation Eq. then becomes 

1 — TjtTk'^ —iark^ iaTk2 
iark^ + St 1 — rjtrk^ —iarki 
~iaTk2 iarki 1 — rjtTk^ 











. BtZ . 







' Bqi ' 




B02 




. Bo3 . 



(62) 
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which is the evolution equation for the a^fi dynamo [l|. We seek solutions of the eigenvalue problem when ^2 = 0. 
Of the three eigenvalues, CTi = (1 — ijtTk'^) is irrelevant, because the corresponding eigenvector does not satisfy the 
solenoidality condition (V • B) = . The remaining two eigenvalues arc 

a± ^ 1 - TTjtk^ ± T{a^k^ -iaksSf^^ , (63) 

corresponding to the eigenvectors 

(^aks , ±i\/ a^k"^ — iak^S , —aki^ . (64) 
The growth rates X± of these eigenmodes are given by 

X± = — ln(fT±) = —ritk^ ± (^a^k"^ — iak^S^^^'^ (65) 
which arc the same as one would get in the case of the a^fl dynamo [l[. 



V. DISCUSSION AND CONCLUSIONS 



This paper presents studies of dynamo action in turbulent shear flows when the turbulence has a non zero correlation 
time. Our goal is to study the dynamics of a system which is complex enough to be a useful model, yet tractable 
analytically; the renovating flows discussed earlier by several authors [l3| (and references therein) provide just such 
a platform. Our contribution is to consider random, helical renovating flows in the context of a background linear 
shear flow. 

We began with a review of the work of Gilbert and Bayly (GB) ^ on random helical renovating flows in the 
absence of a background shear flow. GB considered random flows, each of whose realizations was a plane, sinusoidal 
helical wave. The merit of choosing such simple random ensembles is that the trajectories of fluid elements, in the 
flow caused by each wave, can be integrated analytically. Thus the Green's function mapping the magnetic field from 
one time step to another can be obtained, and averaged over the underlying PDF of the random ensemble of flows. 
GB give the final result, while skipping almost all the intermediate steps. We found it useful not only to record these 
missing steps, but lay them out for the reader so that it becomes easier for us to present our analysis of the more 
complicated problem of renovating flows with shear. 

We then formulated the problem of renovating flows in the presence of a background linear shear flow. Following 
GB, we considered an ensemble of random flows, each of whose members is a plane, sinusoidal helical wave. However, 
unlike in the case considered by GB, the wave cannot be time-independent. In fact, each of these members must be a 
shearing wave, one whose amplitude and wavevector are both time-dependent. Then the trajectories of fluid elements 
(in the flow caused by each shearing wave) were determined analytically, the Green's function (mapping the magnetic 
field from one time step to another) derived, and averaged over the underlying PDF of the random ensemble of flows, 
to obtain a general expression for the (averaged) response tensor. We showed that even without fully averaging the 
response tensor, for which one requires the explicit form of the time-dependent, shearing wave amplitude, it is still 
possible to prove a general result: that there is no dynamo action when the shearing waves are strictly non helical. 

We then considered a particular model in which the shearing waves were generated through external forcing of 
the linear shear flow. Working in the overdamped limit, we derived an explicit form for the response tensor. It is 
interesting to note that this form is closely related to the response tensor of GB; specifically, our response tensor is the 
product of the response tensor of GB with the shearing matrix. This was then applied to the case of a^Q dynamos in 
the limit of small correlation times, and we recovered the standard dispersion relation for the a^n dynamo. Thus the 
growth of the mean field in sheared helical turbulence is as expected from quasilinear closures usually employed to 
derive the mean field equations. This obtains in spite of the fact that magnetic fluctuations in such renovating flows 
are expected to grow more rapidly (as shown explicitly by GB) than the mean field, wherein, one may imagine, that 
quasilinear closures break down. Our work therefore provides another illustration that rapidly growing fluctuations 
need not destroy the growth of the mean field. 

Our result that there is no mean field dynamo in strictly non helical turbulence, even in the presence of shear, raises 
the question as to what causes such growth in numerical simulations @. One possibility is the incoherent a-shear 
dynamo 11 1. There seems to be evidence from such simulations for fluctuations in a (the first reference in Q). Here 
the mean fields are defined as averages over two spatial directions, and fluctuations in a over time are considered as 
meaningful. In the present context, we can define the required fluctuations as due to fluctuations of the parameter 
h from one renovation time to another, which is easier to justify as being physically meaningful. GB themselves 
considered fluctuations in h, but argued that the PDF of h needs to be skewed for there to be net growth. 
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Even in the presence of shear, if we averaged the response tensor say in Eq. (|59p over a mirror-symmetric PDF of 
h with zero mean, the 'hehcal' term would vanish and one would not have a dynamo. However if we think of the 
mean field as being defined before averaging over h, then one could study its dynamics under such fluctuations. In 
the presence of shear, we can see from Eq. (|65p . that both signs of h would cause growth, but the eigenvector in 
Eq. (|64p would get an extra phase shift. It would be interesting to work out exactly how such random changes to the 
eigenvector alters the efficiency of the dynamo. One would expect that a coherent h would lead to a more efficient 
dynamo rather than a fluctuating h. 

The work here has focused on the mean fleld evolution. The same model can also be used to generalize the Kazantsev 
model for the fluctuation dynamo, to the case with finite correlation time and including shear. 
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Appendix A: Calculating the averaged response tensor in the absence of shear 
Below we give the details involved in going from Eq. ([T9|) to Eq. ([2T|) . 



Gij(k) = [5ij + a^gj-r cos(q ■ y + ^j) - hqjhr sm{q ■ y + ip)] e-»k-(aTsin(q-y+^)+br/icos(q-y+^)) ^^1) 

Let the angle between k and q be . Let the component of k perpendicular to q make an angle (f> with a. Now k 
can be written as 



k = 



\^ q(kq) 



g2 



(A2) 



perpendicular to q along q 

Since a and b both lie in the plane perpendicular to q, and also since a and b are perpendicular to each other, 

k-a = k^ • a = |k^||a| cos(^ = (ksin(?)acos(^ 

k • b = ka sin 6 sin (A3) 
Now we average Gij (k) term by term. For the first term, we have 



Ii = e~*'''('^'^*'™('J'y+'/')+'=''^'"=°^(i'y+'^)) (A4) 
The argument of the exponential is written as, 

— irka sin 6 [cos ^ sin(q • y + V') + ^ sin (j) cos(q • y + "0)] = — irkax sin 9 sin(q • y + ■(/; + a) ( A-5) 
where x = (cos^ (j) + h'^ sin^ <PY^^, xcosa = cos0 and xsina = ft, sine/). Hence, we have 



Ii ~ (5jje~*'^'^°x''™^**™(i'y+'''+") (A6) 



First we average over tjj. Since ip goes over all the possible phases, on averaging we get 



2tt 







h = h I e-»(-fcax sine) sin ^ ^ Jo(Tafcxsin0) (A7) 



where we have used the following integral representation of Bessel function of the first kind. 

e^P<^°^^cos{nx)dx = i"7rJ„(/3) (A8) 
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Wc next average over the direction of q. We keep k fixed (say, along z direction) and vary q about k over all the solid 
angles. 



h = Oij / / Jo[TaKxsinff) = / jQ[Takxsmff) smOaO 

Jo Jo 47r 2 7o 

= S,, f Mrakx cos 0) cos 0d0 ^ J„ sin(rafcx) .^gx 
Jo ™fcx 

Lastly we average over the direction of a. Since a can point in any direction in the plane perpendicular to q, we 
average over the angle (j). We get 



sinf 

Ii = Sij goirak, h) , where go{s,h) = (AlO) 

sx 

where the overhead bar now denotes ensemble average over the random variable (jj- We proceed to the remaining term 
of Gy (k) in Eq.dJT]): 



h = [oigj-r cos(q ■ y + V-") - kqjhTsm{q ■ y + t/j)] e-»k-(aTsin(q-y+^)+bT?icos(q-y+^)) (^11) 

We proceed in a similar way as above. We define slightly different parameters. The argument of the exponential is 
now written as, 

— irka sin [cos cj) sin(q ■ y + ^) + h sin <j> cos(q • y + t/j)] = irkax sin cos(q ■ y + — a) (A12) 
where x = {cos"^ cj) + sin^ (f)^^^, xsina = — coscf) and xcosa ^hsin(j) . Therefore we have 



h = [ai^jT cos(q ■ y + V) - hqjhTsm{q ■ y + tp)] e»^fe'^xsinecos(q-y+v-a) (-^13) 
First we average over tp. Since ip goes over all the possible phases, we write 



h = r e'(^'^°>^'^'"'')^°^^ [a.<7jTC0s(C + a) - 6,gj/irsin(C + a)] ^ 

ZTT 



Expanding cos(C + a) and sin(C + a), and keeping only the the even terms under integration, we get 



(A14) 



/ e'(Tfcax sine) cos c [a^g^-r COS a COS C - 6,gj/ir sinacosC] — 
Jo 27r 

IT [a.iqj cos a — bigj h sin a] Ji (jakx sin 0) 

[a,(k • b) - 5,(k • a)] Ji (rafcx sin 6*) (A15) 



ihrqj 



Xak sin 

where we have used Ea. (|A8p to arrive at the second expression and Eq.jASj to obtain the last expression. This can 
be further simplified as 

[ai(k ■ b) - fei(k ■ a)] = [k x (a x b)]^ = ^^[kxa^q]. = qeimnkmqnqj (AI6) 

We next average over the direction of q. We keep k fixed along z direction and vary q about k over all the solid angle 
Then we can write q = sin cos (j) i + sin sin 4> j + cos 0k 

r f^'' ihr . ^ , , ^sin 0d0d4> , 

h = - / — — — a gei„„fc,„g„9j Ji(Tafcxsm6') (A17) 

7o Jo xak sm0 4tt 

The (j) dependence comes only through qn^lj, hence we integrate it first over cj). Note from Eq. ()Al6P that the z 
component of does not contribute to the integral in Eq. (|A17|) . Then we get for the x, y components of qj, 
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Jo 27 



(5n 7 -2 
Hi cm 



sin^e. Thus, 



Sill cue iiLi uqt^^j.^j . 

h — I — ; — : — [a q^iajks] sin 9Ji{Takxs\xv6) ^ — = — / (sin 0Ji(raKx sin 0)) 



g xo-k sin 9 



d 



X2k 

ihraqcij^k^ d 



d{Takx) 
sin{Takx) 



Jq [rakx sin ( 



sin 9d9 ihraqsi^jk^ ^ 
sm9d9 ihraqeij^k^ d 



sin Ode 



X2k d{Takx) Jo 



Jo{Takx sin I 



X2fc d{Takx) 



rakx 



(A18) 



Lastly we average over the direction of a. Since a can point in any direction in the plane perpendicular to q, we 
average over the angle (f). We then get 



l2 = ;7r^^ gi(Tak,h), with gi[s,h) = - [ 



2k 



(A19) 



where the overhead bar now denotes ensemble average over the random variable 
Combining the results for Ii and I2, the response tensor is obtained to be 



GijCk) = Sij go{Tak,h) + 



2fc ri) 



(A20) 



Appendix B: Forced overdamped shearing wave 

The forced, damped Euler equation with a background linear shear is (in this Appendix we use v instead of Uturb 
for brevity): 



d d ^ V 

— + Sxi—]w + Svie2 + (vV)v = -\7p - ^ + f, 



(Bl) 



where tj is a given damping time, and f(x,t) is the external forcing which is assumed to satisfy V • f = 0. The 
pressure, p(x, i), is determined by requiring that Eq. ljBip preserve the incompressibility of the flow. We consider 
external forcing of the form 



f = Re {F(t) exp [iQ(i) • x]} ; Q(i) • F(i) = , 
which excites a single plane shearing wave: 



(B2) 



v = Re{W(t) exp[iQ(i) -x]} 

p = Re{P(0 exp[iQ(i) -x]} 
Incompressibility (i.e. V • v = 0) requires that 



(B3) 



Q-W = 0, 

which makes the nonlinear term, (v • V) v, vanish because 



(B4) 



(W • V) exp [±iQ(t) • x] = ± (iQ ■ W) cxp [iQ(t) • x] = . 
When Eas. (|B3p are substituted in Eq. ljBip , we obtain 



(B5) 



dW f dO \ W 

+ x-^ + 5xiQ2 + 5M^ie2 = -iCtP + F, 

dt \ dt J Td 



(B6) 
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Requiring that the x-dependent terms inside the parentheses vanish imphes that Q(i) must be of the form, 



Qi = qi - Stq2 ; Q2 = 92 ; Q3 
where q = (91,927 93) is a constant wavevcctor. Then W(i) satisfies 



93 : 



(B7) 



dW W 

^ + SW^e2 = -iQP - — 
dt Td 



(B8) 



where = Q • Q = (91 — Stq2)^ + 92 + 93 • We consider the overdamped case when \dW/dt\ ^ |W/rd|, so we drop 
the time derivative term on the left side of Eq. (|B8|) . P can now be cUminated by taking the dot product of Eq. (jB8P 
with Q and using Q • W = 0. Then W(t) satisfies 



SWie2 = SWi ( ^] - — + F. 



The solution is 



Wl = TdFl + STd 



W2 = TdFa - STd 



Td 



Q1Q2 



(B9) 



W2 





- STdQiQ2 




Ql + Ql 




- STdQiQ2 




Q2Q3 




- StaQiQ2 



TdFl 



TdFl 



(BIO) 



Using Q • F = 0, it can be verified that Q • W = 0. The above solution for W is valid for quite arbitrary forms of the 
forcing. Now we make a specific choice for F{t): 



TdF2 = G2+ StGi ; TdFg - G3 



(Bll) 



where G = (Gi, G2, G3) is a constant complex vector that is orthogonal to q (i.e. q • G = ). Then the dependence 
of W on the time t is given in explicit form as 



Wl = Gi + STd 



Q1Q2 



W2 



Wo 



- STdQiQ2 
G2 + StGi - STd 



Gi 



G. 



5Td 



Q^-STdQiQ2 
Q2Q3 



Gi 



- STdQiQ2 



Gi 



(B12) 



We now write the velocity field in explicit real form, using arguments familiar from the discussion of the polarization 
of monochromatic plane electromagnetic waves. G is a complex vector. If its square, G^ = G • G, has argument equal 
to 2ip, then we may write, G = E exp [itp] with q • E = , where E is a complex vector whose square, = E • E, 
is a real quantity. We now express E in explicit form as, E = {he — ia) with q • a = and q • c = , where c 
and a are real vectors orthogonal to q, and /i is a real number; we can choose |c| = |a| and ~1 < h < 1 . Since 



2ih c • a] has been chosen to be a real quantity, we must have c • a = 0. In other words, a and c 



are mutually orthogonal vectors lying in the plane perpendicular to q. Then the velocity field of the sheared plane 
wave of Eq. (jB3[) is given by. 
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v(x, t) = A{t, q) sin [Q(t) -x + tp] + C{t, q) cos [Q(t) • x + -0] , where 



^1,3 = 01,3 + STdai 



Ci,3 = hi ci,3 + St^ci 



2 - STdQlQ2 
Ql,3Q2 



A2 = 02 + Stai — STdCLl 



Ql + Ql 



_Q^~STdQlQ2 

where Qi = qi - Stq2 , Q2 = q2 , Qs ^ qs; 



- STdQiQ2 

C2 = h {C2+ StCl - STdCl 



l + Ql 



q a 



g2 - STdQiQ2 
0, c = qxa, — l</i<l 



(B13) 
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